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SPECIFICATION 

TITLE OF THE INVENTION 

Me thod and apparatus for generating check matrix 

TECHNICAL FIELD 

The present invention relates to a method and an apparatus for 
generating a check matrix for a low-density parity-check (LDPC) code that is 
applied as an error correcting code, and more particularly, to a method and an 
apparatus for generating a check matrix capable of searching a definite and 
characteristic-stabilized check matrix for the LDPC code. 

BACKGROUND ART 

A conventional method of generating check matrixes for LDPC codes 
will be explained below. In a conventional LDPC code encoding/decoding 
system, a communication apparatus at a sending side has an encoder and a 
modulator. On the other hand, an apparatus at a receiving side has a 
demodulator and a decoder. Prior to explanation of the conventional method 
of generating check matrixes for LDPC codes, flows of encoding and decoding 
using LDPC codes are explained first. 

At the sending side, the encoder generates a check matrix H according 
to the conventional method described later. Then, a generator matrix G is 
obtained based on the following condition. 

G: k*n matrix (k: information length, n: code word length) 
GH T =0 (T denotes transposition) 
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The encoder then receives a message (m, m 2 ... m k ) of an information 
length k, and generates a code word C using the generator matrix G. 
C=(mi m 2 ••• m k )G 

=(ci c 2 ... c n ) (where H(ci c 2 ... c n ) T =0) 
5 The modulator subjects the generated code word C to a digital 

modulation such as binary phase shift keying (BPSK), quadrature phase shift 
keying (QPSK), and multi-valued quadrature amplitude modulation (QAM), and 

sends the modulated signal. 

At the receiving side on the other hand, the demodulator receives the 
1 o modulated signal via the channel, and subjects it to a digital demodulation such 
as BPSK, QPSK, and multi-valued QAM. Then the decoder performs an 
iterative decoding by "sum-product algorithm" with respect to the demodulated 
result that is LDPC-coded, and outputs an estimated result (corresponding to 

the original mi m 2 ... m k ). 
i 5 A conventional method of generating check matrixes for LDPC codes 

is explained below. As a check matrix for LDPC codes, the following matrix is 
proposed in a Non-patent Literature 1 (see Fig. 16). 

The matrix shown in Fig. 16 is a binary matrix of "1" and "0", in which a 
part of "1" is hatched. Other parts are all "0". In this matrix, the number of 
20 -1-s in one row (expressed as a row weight) is equal to 4, and the number of 
■1-s in one column (expressed as a column weight) is equal to three. All 
columns and rows have respective uniform weights. Thus, it is generally 
called a "regular-LDPC code". In the codes in the Non-patent Literature 1 , the 
matrix is separated into three blocks, for example, and the second and the third 
25 blocks are subjected to random permutation, as shown in Fig. 16. 



Because the random permutation has no certain rule, it is required to 
execute a time-consuming search by a computer to find codes with a better 
characteristic. 

It is proposed in Non-patent Literature 2 a method using Euclidean 
geometry codes as the LDPC codes that exhibit a relatively stable and 
satisfactory characteristic and can definitely generate a matrix without the use 
of the computer search. This method explains the "regular-LDPC code- 
consisting of regular ensembles. 

The Non-patent Literature 2 proposes a method of generating a check 
matrix for LDPC codes using Euclidean geometry codes EG (2, 2 6 ) as a kind of 
finite geometric codes. This method achieves a characteristic that is located 
closely but 1.45 decibels away from the Shannon limit at an error rate of 10" 4 . 
Fig. 17 is a diagram of a configuration of Euclidean geometry codes EG (2, 2 2 ), 
which has a structure of "Regular-LDPC Codes" with row and column weights 

of 4 and 4, respectively. 

Euclidean geometry codes EG (m, 2 s ) have a characteristic defined as 



n=2 2s -1 



follows: 

Code length: 
Redundant bit length: n-k=3 s -1 
Information length: k=2 2s -3 s 
Minimum distance: d min =2 s +1 
Density: r=2 s /(2 2s -1). 

As can be seen from Fig. 17, Euclidean geometry codes have a 
structure with a cyclically sifted location of "1" in each row from an adjacent row. 
This structure can characteristically configure codes easily and definitely. 



The check matrix generating method in the Non-patent Literature 2 
further includes changing row and column weights based on the Euclidean 
geometry codes to extend rows and columns, if necessary. For example, 
when a column weight in EG (2, 2 2 ) is separated into halves, in the Non-patent 
Literature 2, every other one of four weights located in one column is separated 
into two groups. Fig. 18 is a diagram of an exemplary regular separation of 
the column weight from 4 into 2. 

On the other hand, Non-patent Literature 3 has reported that 
"irregular-LDPC codes" have a better characteristic than that of "Regular-LDPC 
) Codes". This is theoretically analyzed in Non-patent Literature 4 and 

Non-patent Literature 5. The "irregular-LDPC codes" represent such LDPC 
codes that have non-uniformity in either or both of row and column weights. 

Particularly, in the Non-patent Literature 5, a "Sum-Product Algorithm" 
for LDPC codes is analyzed on the assumption that a log likelihood ratio (LLR) 
5 between an input and an output at an iterative decoder can be approximated in 
a Gaussian distribution, to derive a satisfactory ensemble of row and column 
weights. 

Non-patent literature 1: R.G. Gallager, "Low-Density Parity Check 
Codes", M.l.T Press, Cambridge, MA, 1963. 
20 Non-patent literature 2: Y. Kou, S. Lin, and M. P. C. Fossorier, "Low 

Density Parity Check Codes Based on Finite Geometries: A Rediscovery", ISIT 
2000, pp. 200, Sorrento, Italy, June 25-30, 2000. 

Non-patent literature 3: M. G. Luby, M. Mitzenmacher, M. A. Shokrollahi, 
and D. A. Spielman, "Improved Low-Density Parity-Check Codes Using 
25 Irregular Graphs and Belief Propagation", Proceedings of 1998 IEEE 
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. International Symposium on Information Theory, pp. 171 , Cambridge, Mass., 

August 16-21, 1998. 

Non-patent literature 4: T. J. Richardson and R. Urbanke, "The 
capacity of low-density parity-check codes under message-passing decoding", 
5 IEEE Trans. Inform. Theory, vol.47,No.2, pp.599-618, Feb. 2001. 

Non-patent literature 5: S.-Y. Chung, T. J. Richardson, and R. Urbanke, 
"Analysis of Sum-Product Decoding of Low-Density Parity-Check Codes Using 
a Gaussian Approximation", IEEE Trans. Inform. Theory, vol.47, No.2, 

pp.657-670, Feb. 2001. 

1 o According to the conventional method of generating check matrixes for 

LDPC codes disclosed in the Non-patent Literature 5, however, the number of 
"1" points in a row (corresponding to a degree distribution of variable nodes 
described later) and the number of "1" points in a column (corresponding to a 
degree distribution of check nodes described later) are both employed as 

1 5 variables to derive the degree distribution of variable nodes and the degree 
distribution of check nodes that can maximize the following equation (1 ) (rate: 
coding rate). In other words, a linear programming is employed to search an 
ensemble that minimizes a Signal to Noise Ratio (SNR). 

JpW (d 

rate = l — 

to). 

Jo 

20 Therefore, a problem arises that a check matrix derived from the 

maximum of the "rate" has a flux and unstable characteristic. In addition, the 
conventional method of generating check matrixes for LDPC codes iteratively 
executes the derivation of the degree distribution of variable nodes and the 
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derivation of the degree distribution of check nodes over certain times. 
Therefore, a problem arises that it takes time to some extent for searching. 

The present invention has been achieved in consideration of the above 
problems, and accordingly has an object to provide a method of generating 
5 check matrixes for LDPC codes having satisfactory performance capable of 
easily searching a definite and characteristic stabilized check matrix for LDPC 
codes corresponding to an optional ensemble, and a check matrix generating 
apparatus. 

10 DISCLOSURE OF THE INVENTION 

In order to solve the above problems and to achieve the object, a 
method of generating a check matrix according to the present invention 
includes determining a coding rate, generating a basic matrix that satisfies 
conditions that "weights of rows and columns are constant" and "number of 
15 cycles is equal to or more than six", determining number of columns and 

number of rows of the check matrix to be finally obtained, substituting rows of 
the basic matrix created, based on a specific relational equation, searching 
provisionally an ensemble of row weights and column weights of the 
low-density parity check code by executing a Gaussian approximation based 
20 on a predetermined condition before a row deletion, deleting rows of the basic 
matrix after the substituting in order from a bottom by considering the number 
of rows after a division, that is, the number of rows of the check matrix to be 
finally obtained, searching an optimal ensemble of row weights and column 
weights of the low-density parity check code by executing the Gaussian 
25 approximation based on a predetermined condition after the row deletion, and 
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dividing at random the row weights and the column weights of the basic matrix 
after the row deletion based on the optimal ensemble. 

According to the present embodiment, first, a coding rate is determined. 
Next, for example, a "basic matrix based on an integer lattice structure" having 
5 constant weights of rows and columns and a minimum number of cycles as 8 is 
created. The created basic matrix is substituted such that the weights in the 
matrix are arranged at the top of the columns. An ensemble of 
irregular-LDPC codes is provisionally searched by Gaussian approximation 
based on a condition before row deletion. Rows are deleted in order from the 
1 0 bottom of the basic matrix after the permutation by considering the number of 
rows after a division. An optimal ensemble of irregular-LDPC codes is 
searched by Gaussian approximation based on a condition after the row 
deletion. Finally, the weight of the basic matrix after the row deletion is 
divided at random according to a predetermined procedure based on this 
15 ensemble. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 is a flowchart of a method of generating check matrixes for 
LDPC codes according to the present invention; Fig. 2 is a diagram of an 
LDPC encoding/decoding system; Fig. 3 is a diagram of a lattice structure 
when 1 (x, y)=mxx + y + 1 , m-5, and k=3; Fig. 4 is a diagram of m classes using a 
slope s line as a unit; Fig. 5 is a diagram of an algorithm to design an 
aggregate that sets 8 as a minimum number of cycles; Fig. 6 is a diagram of a 
search result when the algorithm shown in Fig. 5 is executed using m=5 and 
25 k=3; Fig. 7 is a diagram of one example of a basic matrix; Fig. 8 is a diagram of 
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a search result when the algorithm shown in Fig. 5 is executed using m=353 
and k-10; Fig. 9 is a diagram of an algorithm of permutation (rearrangement) 
of a basic matrix; Fig. 10 is a diagram of one example of a basic matrix after a 
permutation; Fig. 1 1 is a diagram of a final ensemble of a generator function 
Mx) and a generator function o(x); Fig. 12 is a diagram of a dividing procedure 
according to a conventional paper; Fig. 13 is a diagram of a basic random 
sequence C(i) and permutation patterns of the basic random sequence; Fig. 14 
is a diagram of Latin square sequences L iq (i); Fig. 15 is a diagram of a relation 
between an Eb/No and a bit error rate; Fig. 16 is a diagram of a conventional 
check matrix for LDPC codes; Fig. 17 is a diagram of a configuration of 
Euclidean geometry codes EG (2, 2 2 ); and Fig. 18 is a diagram of an 
exemplary column weight regularly separated from 4 into 2. 

BEST MODE FOR CARRYING OUT THE INVENTION 

Exemplary embodiments of a method and an apparatus for generating 
a check matrix according to the present invention will be explained in detail 
below with reference to the accompanying drawings. It should be noted that 
the embodiments are not intended to limit the invention. 

Fig. 1 is a flowchart of a method of generating check matrixes for 
LDPC codes according to the present invention. The method of generating 
check matrixes for LDPC codes according to the present embodiment may be 
executed within a communication apparatus according to set parameters, or 
may be executed within other control unit (such as a computer) outside of the 
communication apparatus. When the method of generating check matrixes 
5 for LDPC codes according to the present embodiment is executed at the 
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outside of the communication apparatus, generated check matrixes for LDPC 
codes are stored within the communication apparatus. In the following 
embodiment, the above method is executed within the communication 
apparatus for the sake of explanation. 
5 Prior to explanation of the method of generating check matrixes for 

LDPC codes according to the present embodiment, the positions of an encoder 
and a decoder capable of achieving the method are explained first together 
with the conventional method of generating check matrixes for "irregular-LDPC 
codes". 

10 F ig. 2 is a diagram of an LDPC encoding/decoding system. In Fig. 2, 

a communication apparatus at a sending side includes an encoder 101 and a 
modulator 102, and a communication apparatus at a receiving side includes a 
demodulator 104 and a decoder 105. Flows of encoding and decoding using 
LDPC codes will be explained below. 

At the sending side, the encoder 101 generates a check matrix H using 
the method of generating check matrixes for LDPC codes according to the 
present embodiment described later. Then, a generator matrix G is derived 
from the following condition. 

G: k*n matrix (k: information length, n: code word length) 
20 GH T =0 (T denotes transposition) 

The encoder 101 then receives a message (im m 2 ... m k ) of an 
information length k, and generates a code word C using the generator matrix 
G. 

C=(mi m 2 ... m k )G 
25 =(ci c 2 ... c n ) (where H(ci c 2 ... c n ) T =0) 
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The modulator 102 subjects the generated code word C to a digital 
modulation such as BPSK, QPSK, and multi-valued QAM, and sends the 
modulated signal. 

At the receiving side, on the other hand, the demodulator 104 receives 
the modulated signal via the channel 103, and subjects it to a digital 
demodulation such as BPSK, QPSK, and multi-valued QAM. The decoder 
105 then subjects the LDPC-coded, demodulated result to an iterative 
decoding by "sum-product algorithm" and provides an estimated result 
(corresponding to the original mi m 2 ... nn k ). 

The conventional method of generating check matrixes for 
"irregular-LDPC codes" theoretically analyzed in the Non-patent Literature 5 is 
explained next in detail. In this case, a "sum-product algorithm" for LDPC 
codes is analyzed, on the assumption that a log likelihood ratio (LLR) between 
an input and an output at an iterative decoder can be approximated in a 
Gaussian distribution, to derive a satisfactory ensemble of row and column 
weights. 

The method of generating check matrixes for LDPC Codes described 
in the Non-patent Literature 5, or Gaussian Approximation, has a premise that 
defines a point of "1" on a row as a variable node and a point of "1" on a column 
as a check node in the check matrix. 

LLR message propagation from a check node to a variable node is 
analyzed first. The following function (2) is defined on condition that 0<s<co 
and 0<t<oo. In this case, s=m u0 denotes a mean of uO; uO an LLR associated 
with a signal received via a channel containing a Gaussian noise of distributed 
value a n 2 ; and t an ensemble average of LLR output values at check nodes at 
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the time of certain iteration. 



f j (s,t) = <t)- 1 
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f(s,t) = £ Pj f i (s,t) 

j=2 



and P i denote ratios of edges belonging to variable nodes and check 
nodes at a degree of i, respectively, d, denotes a degree of maximum variable 
nodes, and d r a degree of maximum check nodes. X(x) and p(x) denote 
generator functions of degree distribution associated with variable nodes and 
check nodes, and can be represented by the equations (3) and (4), 
respectively. (A degree is expressed as the number of "1"s in each row of 
variable nodes or each column of check nodes). 

(3) 



i = 2 



pM=2> xI " 1 



(4) 



i=2 



where, «x> is defined as shown in the following equation (5). 
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if x<0 



(5) 



The equation (2) can be represented equivalent^ by the following equation (6). 

(6) 

t, =f(s,t,_ 1 ) 

where t, denotes an ensemble average of LLR output values on check nodes at 

the time of the l-th iteration. 

A condition required for deriving an SNR limit (threshold) that provides 
an error with a value of 0 includes t,(s)^ (expressed as FT) when l^o. In 
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order to satisfy this condition, it is required to satisfy the following conditional 
equation (7). 

t<f(s,t), allteR + (7) 
LLR message propagation from a variable node to a check node is 

analyzed next. The following equation (8) is defined on condition that 0<s<co 

and 0<r<1 . In this case, r has an initial value r 0 of Ms). 

's + Ci-l^Pj-ttl-Cl-rH 

j= 2 



hj(s,r) = <|> 



(8) 

h(s,r) = 2 X i h '( s,r ) 

i = 2 

The equation (8) can be represented equivalents by the following 
equation (9). 

r, = h(s,r l _i) ^ 
A condition required for deriving an SNR limit (threshold) that provides 

an error with a value of 0 includes r,(s)->0. In order to satisfy this condition, it 

is required to satisfy the following conditional equation (10). 

r>h(s, r), all re(0, tfs)) ( 10 ) 
In the Non-patent Literature 5, optimal degrees are searched for 

variable nodes and check nodes using the above equation in the following 

procedure (Gaussian Approximation). 

(1) On the assumption that a generator function Mx) and a Gaussian noise 
a n are given, a generator function p(x) is used as a variable to search a point 
that maximizes the equation (1) previously described. A constraint condition 
associated with this search includes normalization to p(x)=1 and satisfaction of 
the equation (7). 

(2) On the assumption that a generator function p(x) and Gaussian noise a n 
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are given (as a value resulted from the procedure (1), for example), a 
generator function X(x) is used as a variable to search a point that maximizes 
the equation (1). Aconstraint condition associated with this search includes 
normalization to Mx)=1 and satisfaction of the equation (10). 

(3) In order to find the maximum "rate", the procedures (1) and (2) are 
iteratively executed to search a better ensemble of the generator function X{x) 
and the generator function p(x) with a linear programming. 

(4) Finally, signal power is normalized to 1 based on Gaussian noise a n to 
find an SNR limit (threshold) (see the equation (11)). 

threshold(dB) = -l 0 * logl 0(2 * o 2 n ) ^ 1 1 ^ 

A problem is found in the Non-patent Literature 5 that a check matrix 
derived from the maximum of the "rate (coding rate)" is flux, and the rate fixed 
in design as a specification varies. In addition, in the Non-patent Literature 5, 
the derivation of the degree distribution associated with variable nodes and the 
derivation of the degree distribution associated with check nodes are iteratively 
executes over certain times. Therefore, a problem arises that it takes time to 
some extent for searching. Further, a problem arises that the check matrix 
does not easily apply to an optional ensemble, an optional code length, and an 

optional coding rate. 

In the present embodiment, a method of easily searching in a short 
time a definite and characteristic-stabilized check matrix for "irregular-LDPC 
codes", which corresponds to an optional ensemble, an optional code length, 
and an optional coding rate, is explained (see Fig. 1). Specifically, the check 
matrix for "irregular-LDPC codes" is generated by using a basic matrix 
described later (definition: "row and column weights are constant", and 



14 



"number of cycles is equal to or more than 6"). 

In the method of generating a check matrix for LDPC codes according 
to the present embodiment, a coding rate is determined first (step S1). As an 
example, the coding rate is set to 0.5. 

A basic matrix based on an integer lattice structure is created on the 
assumption of the basic matrix (definition: "row and column weights are 
constant", and "number of cycles is equal to or more than 6") that is necessary 
to obtain a check matrix of "irregular-LDPC codes" (step S2). In the 
coding/decoding using LDPC codes, generally when the number of "cycle 4" 
and "cycle 6" is as small as possible on a binary graph, satisfactory 
characteristics can be obtained. Therefore, LDPC codes having a 
construction of restricting the occurrence of a small number of "cycle 4" and 
"cycle 6" is preferable. Therefore, in the present embodiment, a basic matrix 
using a minimum number of cycles as 8 is created. A procedure of creating a 
basic matrix (minimum number of cycles as 8) based on an integer lattice 
structure will be explained. A basic matrix using a Cayley graph, a basic 
matrix using a Ramanujan graph, or other matrix can be used so long as the 
basic matrix satisfies the above definition. 

The procedure of creating a basic matrix based on an integer lattice 

structure will be explained. 

(1) First, an aggregate of lines that connect between points (a 
combination of points) in an integer lattice structure is designed. For example, 
an aggregate of a lattice structure is set as L-(x, y). In this case, x is an 
integer of 0<x<k-1 , y is an integer of 0<y<m-1 , k is an integer, and m is a prime 
number. I(x, y) is a linear mapping to map the aggregate L to an aggregate V 
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of points. Fig. 3 is a diagram of a lattice structure when 1(x, y)=mxx+y+1 . 
m=5, and k=3. In Fig. 3, a combination of points that satisfies a certain 
condition is called a line (block). For example, a line having a slope s 
(0<s<m-1) consists of a line of l(x, a+sx mod m) having a point (0, a) as a 
starting point, where a is 0<a<m-1 . As a result, m classes having the line of 
the slope s as a unit are created. Fig. 4 is a diagram of m classes having the 

line of the slope s as a unit. 

(2) An aggregate having a minimum number of cycles as 8 is designed 

based on the algorithm shown in Fig. 5. In other words, a number of columns, 

a number of rows, weight of columns, and weight of rows are searched. In Fig. 

5, S denotes an aggregate of s, and B(s) denotes an aggregate of a class 

corresponding to the slope s. 

Fig. 6 is a diagram of a search result when the above algorithm is 
executed using m=5 and k-3. N-|B| denotes a number of columns in the 
basic matrix, M=|V| denotes a number of rows in the basic matrix, d c denotes 
weight of rows in the basic matrix, and d v denotes a number of columns in the 
basic matrix. Fig. 7 is a diagram of an execution result of the above algorithm, 

that is, a basic matrix. 

The basic matrix needs to be designed in a relatively large size in order 
to realize a division of rows and columns described later. When the above 
algorithm is executed using k-10 and m=353, a search result shown in Fig. 8 is 
obtained. When <fc-1 0 and d c -7 as shown in Fig. 8, a divide processing 
described later can be carried out. In the present embodiment, the 
processing at step S2 is executed to make firm a number of columns N (-2471 ) 
5 in the basic matrix based on the integer lattice structure and a number of rows 
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M (=3530) in the basic matrix based on the integer lattice structure. 

Next, a number of columns N' of a check matrix (a check matrix of 
"irregular-LDPC codes") finally obtained is determined (step S3). At this point, 
a number of rows M' of the check matrix of "irregular-LDPC codes" is 
determined as M'=N'* (1-rate). For example, when N-6000 and rate=0.5, M' 
is determined as M'=6000x0.5=3000. 

Next, the basic matrix created above is rearranged according to an 
algorithm shown in Fig. 9 such that 1s are located at higher positions in the 
columns (step S4). Fig. 9 is a diagram of the algorithm of permutation 
(rearrangement) of the basic matrix. Fig. 10 is a diagram of one example of a 
basic matrix after a permutation. The basic matrix when k=3 and m=5 shown 
in Fig. 7 is rearranged according to the algorithm shown in Fig. 9. Based on 
this permutation, columns of large weight can be left as far as possible and a 
variation in the weight within columns can be reduced as far as possible, in 
delete processing described later. 

Next, an ensemble (degree distribution) of "irregular-LDPC codes- 
based on a requested coding rate is provisionally obtained using optimization 
by Gaussian approximation (step S5). In this case, 7i (1-1 , 2, .... max, 
2< Y i<Y2<...<y m ax) denotes a column weight, p. and u+1 (2<p<d c -1) denote row 
weights, d c denotes a weight of a basic matrix, ^_ Yi (0<^D denotes a ratio 
of an edge that belongs to the column weight Yi , p_p and p_(u+1) (0<" P _p, 
p_(n+1)- *1) denote ratios of an edge belonging to the row weights p and u+1 , 
b_n and b_(n+1) denote nonnegative integers, Mx) denotes a generator 
function of a column weight distribution, p(x) denotes a generator function of a 
row weight distribution, n_p and n_(u+1) denote numbers of rows of the row 
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ghts n and n+1 . and n_ Yi denotes a number of columns of the column 
ght Yi . The above X(x) and p(x) are defined by the equation (13). 

p(x) = p_nxx^ 1 +(l-P_^)xx^ (12) 
b_n x n + b_(n + 1)* (n + 1) = dc 



uxb |i 

p _^ j= ____ (13) 

x (U + I)xb_(n + 1) 

PJP + W- A 

The execution procedure of Gaussian approximation according to the 
present embodiment to search an ensemble of the generator function X(x) of a 
column weight distribution and the generator function p(x) of a row weight 
distribution will be explained. 

(1) A coding rate "rate" is fixed (step S1). 

(2) A generator function A.(x) and a generator function p(x) are 
simultaneously used as variables, and a linear programming is employed to 
search optimal generator functions Mx) and o(x) that can maximize Gaussian 
noise a n (see the following equation (14)). A constraint condition associated 
with this search is to satisfy the equation (18) described later. 



h,(s,r) = <t> 



(•max / ,. ,\ 

(i-i^pMi-Ci-^ ) 



s + l. 

w (14) 



15 h(s,r)=£x i h j (s,r) 

i=2 

!_ 1 Ltanhii-e-iH^ldn, ifx>0 



ifx<0J 



The above s denotes an average of a log likelihood ratio (LLR) 
between a binary signal of {-1, 1} received as a transmission signal and a 
signal received through the Gaussian channel, and can be derived from 
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S=2/a n 2 - 

As explained above, according to the present embodiment, the 
generator functions X(x) and p(x) that satisfy a predetermined condition are 
obtained at one-time linear programming. Therefore, it is possible to create a 
definite and characteristic-stabilized ensemble more easily in a shorter time 
than it is by the method described in the Non-patent Literature 5 that iteratively 
executes derivation of the generator functions X(x) and p(x) to derive both 
optimal values. 

The number of rows M' (=3000) is obtained after executing the divide 
processing of rows described later following b_u, b_(n+1). * and ^+1 obtained 
at step S5. The number of rows M (=3530) of the basic matrix is obtained. 
In this case, rows of a number shown in the following equation (15) are deleted 
in order from the bottom of the basic matrix after the permutation (step S6). In 
this example, (3530-3000)/(1 +0)=530 rows are deleted. A matrix after 
deleting the rows has a set of row weights as {di, d 2 , .... d v >. 

M-M' 1 (15) 
b_|i + b_(n + l) 

Next, an ensemble (degree distribution) of "irregular-LDPC codes- 
based on a requested coding rate is obtained using optimization by Gaussian 
approximation, using the following constraint equations (16), (17), (18), and 

(19) (step S7). A matrix p={2, 3 d v } expressed by (p>j of the equation (16) 

represents a matrix of a nonnegative integer l(p)*L that includes all elements 
satisfying the equation (16). A matrix shown expressed by ((3),,j of the 
equation (17) represents a square matrix of a nonnegative integer l*l. 
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Fig. 11 is a diagram of a final ensemble of the generator functions m 
and p(x) obtained at step S7 when the ensemble is adjusted in the above 
procedure. 

Finally, a dividing procedure of one row or one column in the basic 
10 matrix after the permutation (step 88) will be explained. The Non-patent 
Literature 2 proposes a regular dividing method concerning the dividing 
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procedure. Fig. 12 is a diagram of the dividing procedure in the literature. 
First, a matrix is numbered as shown in Fig. 12. In this example, column 
numbers are given as 1 , 2, 3, and so on in order from the left end, and row 
numbers are given as 1 , 2, 3, and so on in order from the top. For dividing 32 
points*one column into 8 points*4 columns, for example, this is regularly 
divided according to the following equation (20). 

(20) 

S m (n) = B,(m+4*n) 

In the above equation, m=1 , 2, 3, and 4, and n-0, 1,2,3, 4, 5, 6, and 7, 
and I denotes a column number of EG (2,2 5 ). B,(x) denotes a position of "1" in 
an l-th column of EG (2,2 5 ). S m (n) denotes a position of "1" in an m-th column 

of the matrix after the division. 

Specifically, a row number that shows a position of "1" in the l-th 
column of EG (2,2 5 ) is as follows. B l( x H 1 32 114 136 149 223 260 382 402 
438 467 507 574 579 588 622 634 637 638 676 717 728 790 851 861 879 947 
954 971 977 979 998}. As a result, a row number that shows a position of "1" 
in the first to fourth columns of the matrix after the division is as follows, based 
on a regular extraction of the number of "1" from Bi(x). 
Si(n)={1 149 402 574 634 717 861 971} 
S 2 (n)={32 223 438 579 637 728 879 977} 
S 3 (n)={114 260 467 588 638 790 947 979} 
S 4 (n)={136 382 507 622 676 851 954 998} 
In other words, 32 pointsxone column is divided into 8 points*4 columns. 

On the other hand, in the divide processing of the basic matrix after the 
permutation according to the present invention, a regular dividing like the 
S above processing is not executed but the number of "1" is extracted at random 
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from Bi(x) (see a detailed example of a random division described later). Any 
method of extraction processing can be used so long as randomness is 
maintained. 

When a position of "1" in the m-th column of the matrix after the 
division is R m (n), for example, R m (n) becomes as follows. 
Rl (n)=(1 114 574 637 851 879 977 979} 
R 2 (n)={32 136 402 467 588 728 861 971} 
R 3 (n)={149 260 382 438 579 638 717 998} 
R4(n)={223 507 622 634 676 790 947 954} 

One example of a random division, that is, "a division method using a 
Latin square of random sequences", will be explained in detail. A random 
sequence to carry out the random division is created easily and definitely. An 
advantage of this method is that the same random sequence can be created at 
the sending side and the receiving side. This is extremely important in a real 
system. This method also has an advantage that a condition of code 
characteristic can be accurately prescribed. 

(1) Creation of basic random sequences: 

An example of random sequence creation is described below using 
Euclidean geometry codes EG (2,2 5 ) for convenience of explanation. In 
Euclidean geometry codes EG (2,2 5 ), the number of "1"s present in a row is 
equal to 2 5 =32. 

When P is used for the minimum prime number that satisfies P>d v =2 s , 
for example, P=37 in the case of d v =2 5 . A basic random sequence C(i) with a 
sequence length, P-5=32, is created in accordance with the equation (21). d v 
denotes a maximum weight of a column. Therefore, when codes other than 
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Euclidean geometry codes are selected as a basic matrix, the use of d v of this 
basic matrix makes it possible to apply this division. 
C(1)=1 

C(i+1)=G 0 xC(i) mod P ( 21 > 
5 where, i=0, 1 , .... P-2; and G 0 denotes an original source of Galois Field GF(P). 
As a result, C(i) is represented by the following equation: 

C (i)={1 2 4 8 16 32 27 17 34 31 25 13 26 15 30 23 9 18 36 35 33 29 21 
5 10 20 3 6 12 24 11 22 7 14 28 19}. 

(2) Numbers larger than 32 are deleted so as to obtain the sequence 

10 length d v =2 5 =32. 

C(i)={1 2 4 8 16 32 27 17 31 25 13 26 15 30 23 9 18 29 21 5 10 20 3 6 

12 24 11 22 7 14 28 19} 

(3) A permutation pattern LBj(i) is created using the following equation 

(22). 

15 j-1, 2, •■, d v 

M.2.-.FM < 22 > 
A number larger than LBj(i) is deleted. Fig. 13 is a diagram of the basic 
random sequence C(i) and the permutation pattern LB,(i) of the basic random 
sequence. 

20 (4) A j-th Latin square matrix L ip (i) in the column q and the row i is 

calculated from the following equation (23) to execute a division. When the 
column weight d_p is d_p<d v , numbers larger than d_p are thinned from the 
elements of L iq (i), based on the delete processing at step S6. 
L ip (i)=C(LBj(((q+i-2) mod d v )+1) 

25 2, -, d v 
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i=1,2, -, d v 

q=1,2,-,P-1 < 23 > 
Fig. 14 is a diagram of Latin square sequences L iq (i). The Latin 
square sequences L jp (i) are used to determine a division pattern of the 
5 jx32+q-th column of the matrix to be expanded. For example, the 670-th 
column g 67 o (1) of EG (2,2 5 ) to be shortened by deletion is determined as 
follows: 

g 670 (1)={28 48 84 113 153 220 225 234 268 280 283 284 322 363 374 
436 497 507 525 593 600 617 623 625 644 670 701 783 805 818 892 929}. 
10 This is divided into five columns of the weight 6 and one column of the weight 2. 
The corresponding Latin square L iq (i) is 20x32+30=670. Therefore, the 
following Latin square is obtained: 

L21.30 (1H13 19 9 10 16 24 25 28 23 5 8 12 31 14 30 21 4 6 17 7 15 29 

2 3 27 22 26 18 1 20 32 11}. 
15 As a result, the division pattern becomes as follows: 

g670,l(1) = 9670(L 2 1,30 (1)) 

={322 525 268 280 436 625} i=1 , 2, ■ ■-, 6 

g670,2(1) = g670(L21,30 (1)) 

={644 783 623 153 234 284} i=7,8, ■-, 12 

20 g670,3(1) = g670(L 2 1.30 (1)) 

={892 363 818 600 113 220} 1=13, 14, ■ ■, 16 

g670,4(1)=g670(L 2 1,30 (1)) 

={497 225 374 805 48 84} i=17, 18, -, 24 

g670,5(1) = 9670(L21,30 (1)) 

25 ={701 617 670 507 28 593} i=25, 26, -, 30 
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g670,6(1) = g670(l-21,30 (1)) 

={929 283} i=31,32. 
In general, the elements L j>q (i) of Latin square for g c . e (l) are determined based 
on the following equation (24). 

j = rc/d v l (24) 

g = ((c-l) mod d v )+l 

Characteristics of the above LDPC codes will be compared below. 
Fig. 15 is a diagram of a relation between an Eb/No (a signal power to noise 
power ratio per one information bit) and a bit error rate (BER). A decoding 
method is "Sum-Product algorithm". This characteristic uses the ensemble 
shown in Fig. 11. Fig. 15 is a characteristic comparison between the 
execution of the regular division as described in the Non-patent Literature 2 
and the execution of the division according to a Latin square of random 
sequences. 

As is clear from Fig. 15, according to the regular division as described 
in the Non-patent Literature 2, a large improvement cannot be expected even 
with "irregular-LDPC codes". In contrast, the random division of the present 
embodiment can provide a remarkably improved performance when it is 
implemented because the probability of the occurrence of a loop decreases 
substantially. 

As explained above, according to the present embodiment, first, a 
coding rate is determined. Next, a "basic matrix based on an integer lattice 
structure" having constant weights of rows and columns and a minimum 
number of cycles as 8 is created. The created basic matrix is substituted 
based on a specific relational equation. An ensemble of irregular-LDPC 
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codes is provisionally searched by Gaussian approximation based on a 
condition before row deletion. Rows are deleted in order from the bottom of 
the basic matrix after the permutation by considering the number of rows after 
a division. An optimal ensemble of irregular-LDPC codes is searched by 
5 Gaussian approximation based on a condition after the row deletion. Finally, 
the weight of the basic matrix after the row deletion is divided at random 
according to a predetermined procedure based on this ensemble. With this 
arrangement, a definite and characteristic-stabilized check matrix for 
"irregular-LDPC codes" can be generated easily in a short time corresponding 
10 to an optional ensemble, an optional code length, and an optional coding rate. 

As explained above, according to the present invention, first, a coding 
rate is determined. Next, a basic matrix having constant weights of rows and 
columns and a number of cycles equal to or larger than 6 is created. The 
created basic matrix is substituted based on a specific relational equation. An 
1 5 ensemble of irregular-LDPC codes is provisionally searched by Gaussian 

approximation based on a condition before row deletion. Rows are deleted in 
order from the bottom of the basic matrix after the permutation by considering 
the number of rows after a division. An optimal ensemble of irregular-LDPC 
codes is searched by Gaussian approximation based on a condition after the 
20 row deletion. Finally, the weight of the basic matrix after the row deletion is 
divided at random according to a predetermined procedure based on this 
optimal ensemble. With this arrangement, there is an effect that a definite 
and characteristic-stabilized check matrix for "irregular-LDPC codes" can be 
generated easily in a short time corresponding to an optional ensemble, an 
25 optional code length, and an optional coding rate. 
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INDUSTRIAL APPLICABILITY 

As explained above, the method of generating check matrixes for 
LDPC codes and the check matrix generating apparatus according to the 
present invention are useful for the communication system that employs the 
LDPC codes as error correcting codes. Particularly, the method and the 
apparatus are suitable for a communication apparatus that searches definite 
and characteristic-stabilized LDPC codes. 



